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Abstract 

We present a covariant density matrix approach to kinetic theory of QED plasmas 
subjected to a strong external electromagnetic field. A canonical quantization of 
the system on space-like hyperplanes in Minkowski space and a covariant general- 
ization of the Coulomb gauge is used. The condensate mode associated with the 
mean electromagnetic field is separated from the photon degrees of freedom by a 
time-dependent unitary transformation of both, the dynamical variables and the 
nonequilibrium statistical operator. Therefore even in the case of strong external 
fields a perturbative expansion in orders of the fine structure constant for the cor- 
relation functions as well as the statistical operator is applicable. A general scheme 
for deriving kinetic equations in the hyperplane formalism is presented. 
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1 Introduction 



In recent years the theoretical study of dense relativistic plasmas is of in- 
creasing interest. Such plasmas are not only limited to astrophysics, but can 
nowadays be produced by high- intense short-pulse lasers [1-3]. In view of the 
inertial confinement fusion, one has to consider a plasma under extreme con- 
ditions which is created by a strong external field. This new experimental 
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progress needs a systematic approach based on quantum electrodynamics and 
methods of nonequihbrium statistical mechanics. 

Considerable attention has been focussed on a mean- field (Vlasov-type) kinetic 
equation for the fermionic Wigner function, which is an essential step towards 
transport theory of laser-induced QED plasmas. Using the Wigner operator 
defined in four-dimensional momentum space [4-6], a manifestly covariant 
mean-field kinetic equation can be derived from the Heisenberg equations of 
motion for the field operators. In this approach, however, it is difficult to 
formulate an initial value problem for the kinetic equation since the four- 
dimensional Fourier transformation in the covariant Wigner function includes 
integration of two-point correlation functions over time. This difficulty does 
not appear in the scheme based on the one-time fermionic Wigner function 
where the field operators are taken at the same time and only the spatial 
Fourier transformation is performed. In the context of QED, the one-time 
formulation was proposed by Bialynicki-Birula et al. [7] (referred to in the 
following as BGR) and used successfully in their study of the electron-positron 
vacuum. Within this approach one can explore a number of attractive features. 
The one-time Wigner function has a direct physical interpretation and allows 
to calculate local observables, such as the charge density and the current 
density. The description in terms of one-time quantities is quite natural in 
kinetic theory based on the von Neumann equation for the statistical operator 
and provides a consistent account of causality in collision integrals. 

It should be noted, however, that the one-time Wigner function does not 
contain complete information about one-particle dynamics; the spectral prop- 
erties of correlation functions can be described only in terms of two-point 
Green's functions which are closely related to the covariant Wigner function. 
Recently this aspect of relativistic kinetic theory was studied within the mean- 
field approximation [8,9]. The aforementioned incompleteness of the one-time 
description is well known in non-relativistic kinetic theory, where two-time cor- 
relation functions can, in principle, be reconstructed from the one-time Wigner 
function by solving integral equations which follow from the Dyson equation 
for nonequihbrium Green's functions [10]. The reconstruction problem in rel- 
ativistic kinetic theory remains to be explored. The sohition of this problem 
requires a further development of the relativistic density matrix method as 
well as the relativistic Green's function technique. 

In this and subsequent works we develop a density matrix approach to kinetic 
theory of QED plasmas subjected to a strong electromagnetic field. From the 
conceptual point of view, our aim is to generalize the BGR scheme [7] in two 
aspects. First, we wish to present the one-time formalism in covariant form. 
This removes a drawback of the BGR theory which is not manifestly covariant. 
Second, we will develop a scheme which allows to go beyond the mean-field 
approximation, including dissipative processes in QED plasmas and the inter- 
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play between collisions and the mean- field effects. Whereas subsequent studies 
will be concerned with explicit liinetic equations, the present first part consid- 
ers some general problems of the one-time covariant approach to relativistic 
kinetic theory. In comparison to QED where the main object is the (S'-matrix 
constructed from vacuum averages of the field operators, kinetic theory of 
QED deals with averages over a nonequilibrium ensemble describing a many- 
body system. Therefore we use the Hamiltonian formalism which is typical 
for the density matrix method. In this case, however, one meets with some 
fundamental problems which are considered in this paper. In order to make 
the theory manifestly covariant, canonical quantization of the system will be 
carried out in a covariant fashion. Another point is that in the presence of 
a strong electromagnetic field, perturbation expansions in the fine structure 
constant are not suitable. To overcome this difficulty, we will present a proce- 
dure which allows to separate the classical part of the electromagnetic (EM) 
field and the photon degrees of freedom at any time. 

The paper is organized as follows. In Section 2 we briefiy sketch a scheme of 
relativistic statistical mechanics in the form adapted to kinetic theory. In our 
approach we use a manifest covariant Schrodinger picture on space-like hyper- 
planes in Minkowski space. Analogous formulations of relativistic quantum 
mechanics and quantum field theory can be found in literature for various ap- 
plications (see, e.g., [11-15]). In this way, "equal-time" correlation functions 
are defined with respect to the "invariant time" variable on a hyperplane. 
In Section 3 we perform canonical quantization of QED on space-like hyper- 
planes and derive the covariant quantum Hamiltonian. Section 4 deals with 
the condensate mode which corresponds to the electromagnetic field induced 
by the polarization in the system. The condensate mode is eliminated by 
a time-dependent unitary transformation of the statistical operator and dy- 
namical variables. As a result, we obtain the effective Hamiltonian, where 
the interaction of fermions with the mean electromagnetic field is incorpo- 
rated non-perturbatively at any time, while the interaction between fermions 
and photons is described by a term which can be taken into account within 
perturbation theory. It is shown how Maxwell equations for the mean electro- 
magnetic field are recovered in our scheme. In Section 5 the covariant one-time 
Wigner function and the photon density matrix are introduced and a method 
for deriving kinetic equations in the hyperplane formalism is outlined. The 
paper is summerized with a short discussion of the results and an outlook to 
further applications. 

We use the system of units with c — h — 1. The signature of the metric tensor 
is (+,-,-,-). 
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2 Nonequilibrium statistical operator in the hyperplane formalism 

2.1 The relativistic von Neumann equation 

It is well known that in the special theory of relativity a quantum state of a 
system is defined by a complete set of commuting observables which can be 
associated with a three-parameter space-like surface a in Minkowski space. 
Among these surfaces three-dimensional hyperplanes are especially simple to 
deal with [11-14]. Since the use of arbitrary space-like surfaces does not lead 
to new physics, in what follows we restrict our consideration to hyperplanes. 
A space-like hyperplane a = cr^ is characterized by a unit time-like normal 
vector n^^ and a scalar parameter r which may be interpreted as an " invariant 
time" . The equation of the hyperplane reads 

X ■ n = T, = n^n^ = 1. (2.1) 

In the special Lorentz frame where n'^ = (1, 0, 0, 0) and consequently Eq. (2.1) 
reads — t the parameter r coincides with the time variable t — x°. We will 
refer to this special frame as the "instant frame" , since only here observables 
are measured at the same instant of time t. By treating a state vector [^[c^^^]) 
as a functional of ^, the covariant Schrodinger equation can be derived from 
the relation between the state vector on the hyperplane a and the state vector 
on the hyperplane a' — La which is obtained by an inhomogeneous Lorentz 
transformation L — {a, A}: 

a ^ a' = La : x ^ x' = Ax + a. (2.2) 

The relation between the state vectors is [16] 

U{L)mLa]) = ma]), (2.3) 

where U{L) — U{a,A) is a unitary representation of the inhomogeneous 
Lorentz group. The generators of this representation, and M^", are the 
energy-momentum vector and the angular momentum tensor, respectively. For 
our purpose, the only transformations of relevance are pure time-like trans- 
lations which change the value of r. Recalling the form of U (a. A) for pure 
translations 

C/(a,l)=exp{iPX} (2-4) 
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and introducing the notation = |^(n, r)), Eq. (2.3) can be written 

for an infinitesimal time-like translation — 5t as 

|^(n, r + 6t)) + iSr (P^n'') |^(n, r)) = |^(n, r)) , (2.5) 



from which we obtain the relativistic Schrodinger equation 

i^\^{n,T)) = H{n)\^{n,T)) (2.6) 



with the Hamiltonian on the hyperplane given by 

H{n) = PX- (2-7) 

In the presence of a prescribed external field, the energy-momentum vector 
and, consequently, the Hamiltonian H'^{n) can depend explicitly on r. Com- 
bining Eq. (2.6) with the adjoint equation for the bra-vector, one finds that the 
statistical operator g{n, r) for a mixed quantum ensemble obeys the equation 



dg{n,T) _ . 



gin,T),H^{n) =0, (2.8) 



dr 

which is analogous to the non-relativistic von Neumann equation. 

2.2 Schrodinger and Heisenberg pictures on hyperplanes 

The evolution of a mixed ensemble on space-like hyperplanes can be repre- 
sented in different pictures. The statistical operator in the Heisenberg picture 
does not depend on the parameter r and is associated with some fixed hyper- 
plane Dynamical variables are represented by operators O/^ ([o"„ ^]) which 

are functionals of the hyperplanes. Local operators Ojj{x), which depend on 
the space-time point x, are of particular interest in quantum field theory. In 
what follows it will be convenient to treat such operators as functions of the 
parameter r. To define this dependence, we introduce the transverse projector 
with respect to the normal vector n'*, 

A^, = - n^n„ (2.9) 

and notice that a space-time four-vector x^ can be represented in the form 

x" ^ n^T + xi^, T^n-x, (2.10) 
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where 



(2.11) 



is the transverse (space-hke) component of x. GcometricaUy, Eq. (2.10) means 
that the space-hke vector x'i hes on the hypcrplane ^ passing through the 
space-time point x. Using the decomposition (2.10), a local Heisenberg oper- 
ator Ojj{x) can be written as 

0^(a;) = 0^(nr + a;J ^0^(r,a;J. (2.12) 

Let us assume that P'^ does not depend explicitly on r. Then, recalling the 
well-known equation of motion for Heisenberg operators 

d,OH{x)^-i[dH{x),P,], (2.13) 

one readily finds that the time-like evolution of such operators is described by 
the equation 

0^(t, xJ = S-'-^omn) Q^^^^^ ^^i(r-r,)H{n) (2.14) 

with the Hamiltonian (2.7). The generalization of Eq. (2.14) to situations in 
which the Hamiltonian H'^ depends explicitly on r is obvious. Defining the 
evolution operator U{T,T';n) as the ordered exponent 

U{t, t'; n) = 7; exp i^-i J H^{n) dfj , (2.15) 

we have 

Oh{^, x±) = U\t, To; n) OHirQ, xj U (r, Tq; n). (2.16) 

In the Schrodinger picture, the statistical operator g{n, r) is r-dependent and 
its time-like evolution is governed by Eq. (2.8), whereas operators Og are 
defined on a fixed hyperplane. Assuming the Heisenberg and Schrodinger pic- 
tures to coincide on the hyperplane (T,„^ , Eq. (2.16) implies that the transition 
from the Schrodinger picture to the Heisenberg picture is given by 

Onir, xj = C/t(r, r^; n) ds{xj U{r, r^; n). (2.17) 
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The mean values 0{x) of local dynamical variables can be calculated in both 
pictures. Using a formal solution of Eq. (2.8) 

Q{n, t) = C/(r, To; n) Q{n, Tq) U^{t, Tq; n), (2.18) 



we find that 

0{x)^{6^{r,xjyo^{ds{xjy. (2.19) 

Here and in what follows the symbol (• • -y stands for averages calculated with 
the statistical operator g{n,r). In many problems one is dealing with partial 
derivatives d^O{x) which enter the equations of motion for local observables. 
In the hyperplane formalism, it is convenient to express the partial derivatives 
in terms of the derivatives with respect to r and Xj^. RecaUing Eqs. (2.10) 
and (2.11), we write 

d hV, V^A'^a^A'' . (2.20) 

Then, in the Heisenberg picture, Eq. (2.19) yields the equation of motion 

d,0{x)^{d^dH{r,xjy\ (2.21) 



where 



d . 



-in^ 



Onir, xj, Hnin, t)1 + V^6^(t, xJ. (2.22) 



In the Schrodinger picture, the r-dependence of the mean values appear through 
the statistical operator which obeys the von Neumann equation (2.8). In this 
picture the equation of motion, bearing a formal resemblance to Eq. (2.21), is 
obtained from Eq. (2.19) 



d,0{x) = (d^OsixJ) (2.23) 



with the analogous definition of the operator acting on local dynamical 
variables: 



9^,0 six J = 



-m, 



ds{xJ,H{n)]+V^ds{xJ. 



(2.24) 
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2.3 "Equal-time" correlation functions 

Describing the evolution of the system in terms of hyperplanes, we can in- 
troduce "equal-time" correlation functions of local dynamical variables with 
respect to the invariant time r. Let Oifj{x), 02h{x), . . . , Oj^^{x) be some local 
Heisenberg operators. Then the "equal-time" correlation function for these 
operators can be defined as 

Fi-ki^i^^ Xk±; t) = (Oi^^(xi) • • • Of^ffi^kW", (2.25) 

where n ■ x^ = n ■ X2 = ... = n ■ x,^ — t. In the Schrodinger picture this 
correlation function takes the form 

Fi...k{xi^, x^^; n, r) = (OisK^) • • • O^^sKx))^- (2-26) 
The covariant von Neumann equation (2.8) yields the equations 

^^i...fe(^ix, ■ ■ ■ , ^fex; n, t) = -i([Oi5Kx) • • • 0,sM, H^inW (2-27) 

which can serve as a starting point for constructing the quantum hierarchy 
for the "equal-time" correlation functions. 



3 Hamiltonian of QED on hyperplanes 

We will now apply the foregoing scheme to a relativistic system of charged 
fermions interacting through the EM field. For definiteness, we take these 
fermions to be electrons and positrons, so that protons will be treated as a 
positively charged background which ensures electric neutrality of the system. 
There is no difficulty in describing protons by an additional Dirac field. Having 
in mind applications to relativistic plasmas produced by high-intense short- 
pulse lasers, we assume the system to be subjected into a prescribed external 
EM field which is not necessarily weak. 

3.1 The Lagrangian density 

The first step in formulating the kinetic theory of QED plasmas is to construct 
the Hamiltonian H{n). We start with the classical Lagrange density 

C{x) = Cu{x) + Cem{x) + Ant(a;) + Axt(a^)> (3-1) 
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where jCjj{x) and Cj^^j{x) arc the Lagrangian densities of free Dirac and EM 
fields respectively. £i„t(.x') is the interaction Lagrangian density, and the term 
'^ext(^) describes the interaction of fermions with the external electromagnetic 
field. In standard notation (see, e.g., [17]), we have 



Cj,{x) = i>{x) (^-^^d^ - m) ^l;{x), (3.2) 

^em{x) = -\f,.{x)F^''{x), (3.3) 
Ant(^) = -3^>^)A>'i:,), (3.4) 
Axt(^) = -J,{x)A^.A^), (3.5) 

where = — d^- In the following the electromagnetic field tensor is 
taken in the form = 9^A^ — d^,A^. The current density four- vector will 
be expressed as = e'ljj'y'^ijj with e < 0. We wish to remark that in our 
approach the four-potential of the EM field is decomposed into two terms. The 
variables yl^(a;) correspond to the EM field caused by charges and currents in 
the system, while A^^^(x) is a prescribed external field. In what follows, only 
the dynamical field A'^{x) will be quantized. 



3.2 Canonical quantization on hyperplanes 



A canonical quantization implies that some gauge fixing condition is im- 
posed on A^. For many-particle systems studied in statistical mechanics, the 
Coulomb gauge seems to be the most natural. However, the disadvantage of 
this gauge is that it is not manifestly covariant. Therefore we will use a general- 
ization of the Coulomb gauge condition which is consistent with the covariant 
description of evolution in terms of space-like hyperplanes. To formulate this 
condition, we introduce for any four-vector the decomposition into the 
transverse and longitudinal parts by 

V ^n'^V^^ + Vj^, Vj^^A'^^V'', (3.6) 

where A'^^^ is the projector (2.9). Then a natural generalization of the Coulomb 
gauge condition reads 

V„^^ = 0. (3.7) 



In the special frame where — (1, 0, 0, 0) and A^ = (^4°, A), Eq. (3.7) reduces 
to V • A = 0, which is the usual Coulomb gauge condition. 
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To define canonical variables for the electromagnetic field on a hyperplane 
(Jnri fi^s^ perform the decomposition (3.6) of the field variables and the 
decomposition (2.20) of the derivatives in the Euler-Lagrange equations 



dC 



-d. 



dC 



0. 



(3.8) 



A simple algebra shows that these equations are equivalent to 



dC d dC 



OA, dr \ dA„ 



- V. 



dC 
dAl 



d 



dr XdA", 



'^^(V.A,,) 

dC\ ^ dC 
-V, 



0, 



a(v,A:^), 



0, 



(3.9) 
(3.10) 



where we use the notation f = df /dr for derivatives with respect to r. Equa- 
tion (3.9) allows to eliminate the variable A^^ in the Lagrangian. First we 
rewrite expressions (3.3) and (3.4) in terms of A^^ and A'i using the decom- 
position procedure for the derivatives and the field A'^. As a result we obtain 
the Lagrangian density in the form 



where we have introduced the notation 



(3.11) 



(3.12) 



Note that the last two terms in Eq. (3.11) do not contain A^^ and A'i. Now 
using the expression (3.11) to calculate derivatives in Eq. (3.9) and taking into 
account that, according to the gauge condition (3.7), V^^x = 0, we get 

V,V'^^„=J,, (3.13) 



In the "instant frame", where — (1,0,0,0), this reduces to the Poisson 
equation for ^4°. The solution of Eq. (3.13) is 

A(^'^x) = / da'G{x^-x'Jj^^{r,x'J, (3.14) 
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where the Green function G{x^} satisfies the equation 



(3.15) 



with the three-dimensional delta function on a hyperplane cr„ defined as 

S'{xJ=J -^.^-''""^{P-n). (3.16) 



The solution of Eq. (3.15) for G{x_^) is given by 



The variable A^^ can now be eliminated in the Lagrange density (3.11) with 
the aid of Eq. (3.14). Terms like V^(- • •) can be dropped since they do not con- 
tribute to the Lagrangian L = j Cda under appropriate boundary conditions. 
Then a straightforward algebra leads to 



+ + 4xt -\j da'j^^{r,xjG{x^ - x'Jj^^{r,x'J. (3.18) 

We will treat the fields At as dynamical variables for the EM field and follow 
the Dirac version of canonical quantization of theories with constraints [18,19]. 
The gauge condition (3.7) is one of the constraint equation in this scheme. 
Another constraint equation follows directly from the definition of transverse 
four- vectors, Eq. (3.6), and reads 

nA^^(x)=0. (3.19) 



We now define canonical conjugates for the field variables A'i by 



Obviously the II's are not independent variables since they satisfy the con- 
straint equations V^IIx = 0, and n^U'i = 0. Thus, we have four constraints 
imposed on the canonical variables. Following the standard quantization pro- 
cedure [19], the commutation relations for the field operators and 11^^ can 
be derived. As shown in Appendix A, these commutation relations are 
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where 



(fix 



(2^ 



-ip-(x—x') 



d{p ■ n) 



A"" - 



M u 
PlPj_ 

pI 



(3.21) 
(3.22) 

(3.23) 



In Appendix B the anticommutation relations for the Dirac field operators on 
hyperplanes are derived. The result can be written as 



(3.24) 
(3.25) 



where a, a' are the spinor indices. The matrix 7||(n) is introduced through the 
following decomposition of the Dirac matrices 7'*: 

r = n%{n) + YAn), iM) = ^.7^ YAn) = {S^. - n^n,) (3.26) 



In the special Lorentz frame where = {t,r) and = (1,0,0,0), we have 
7ii = 7° and S^{x^^ — x'^) = 5{r — r'), so that Eq. (3.24) reduces to the well- 
known anticommutation relation for the quantized Dirac field. 



3.3 Derivation of the Hamiltonian 



The classical Hamiltonian on the hyperplane cr^ can be derived in two ways. 
Following the canonical procedure, H'^{n) is obtained by the Legendre trans- 
formation 

H{n)^ j dcrjn^^i'^ + 7fV' + 'i^7r-£|, (3.27) 



where C is given by Eq. (3.18). To find explicit expressions for the variables 
TT and 7f, which arc conjugates to the fields ip and ■?/', we rewrite the Dirac 
Lagrangian density (3.2) using the decomposition (2.20) of derivatives: 



D 



(3.28) 
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Then we have 

Substituting expressions (3.28) and (3.29) into Eq. (3.27), we arrive at the 
classical Hamiltonian. Another way is to start from the classical analog of 
Eq. (2.7) which reads 

H{n)^PX= I dan^T^X, (3-30) 

where T^,^(.x') is the energy- momentum tensor. In order that the quantized 
Hamiltonian be hcrmitian, the energy-momentum tensor must be real. For 
instance, one can use the so-called Belinfantc tensor [20]. When applied to 
the Lagrangian (3.1), the standard derivation of the Belinfante tensor (see, 
e-g-, [17]) gives 

+ F^^xK + (7. +7^ ^ - ^ (-^-^M + h^'^) ■ (3-31) 

Separating the longitudinal and transverse components with respect to the 
normal vector and then eliminating the r-derivatives of the fields with the 
aid of Eqs. (3.20) and (3.29), the classical Hamiltonian on the hyperplane is 
obtained from Eq. (3.30). It can be verified that in both cases we have the 
same expression for H'^{n). The final step is to replace the canonical variables 
A'i, Hx, ip, ■0 by the corresponding quantum operators. As a result, we find the 
Hamiltonian in the form 

H-{n) = HM + H^Min) + 4ntH + 4xtH, (3-32) 

where Hjj{n) and H^j^{ri) are the Hamiltonians for free fermions and the 
polarization EM field respectively, H^^^{n) is the interaction term, and Hl^^{n) 
describes the external EM field effects. In the Schrodinger picture the exphcit 
expressions for these terms are 

^jda^ l^-'-rAn) V, +m) V^, (3.33) 

HEMin) = jda (^F,^.Fr - ^n.^H^^) , (3.34) 
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AntW = / daj^^A^l^ ^ 2 / ^^/ ~ x'±)j\\ix'±)^ (3-35) 
4;tW = / daj^{xJAUr,xJ. (3.36) 
In the EM field Hamiltonian (3.34) the tensor operator 

pi^u ^ yM^^^ _ v'^i'f (3.37) 

contains only the transverse components of the field operators which are 
decomposed as 

> = n'^i,! + i'^. (3.38) 

The longitudinal part A^^ has been eliminated in the interaction Hamilto- 
nian (3.35) by the equation 

V^V^i,, = j„ (3.39) 

which is analogous to Eq. (3.13). As usual, in Eqs. (3.33) - (3.36) normal or- 
dering in operators is implied. The self-energy contribution to the last term 
in Eq. (3.35) is omitted, so that the product : j^^^x^^ :: j^^{x'^) : is under- 
stood. For simplicity, we have written the Hamiltonian for the case that the 
fermionic subsystem is described by one Dirac field. The generalization to a 
many-component case is obvious. 



4 The condensate mode of the EM field 

An essential feature of the dynamical evolution of QED plasmas in a strong 
external field is that the mean values of the canonical operators[^. At = (A'i) 
and Hx = (Hx), just as the mean values of creation and annihilation bosonic 
operators a and related to the canonical operators by plane-wave expan- 
sions, are not zero. Furthermore, they are macroscopic quantities associated 
with the mean EM field induced by the polarization in the system. In the 
language of statistical mechanics, the variables A^^lx) and Il±{x) describe a 
macroscopic condensate mode of the EM field. This fact does not allow to 
apply perturbation theory directly to the Hamiltonian (3.32) because the in- 
teraction of fermions with the condensate mode or, what is the same, with the 

Prom this point onwards symbols A't, 11'^, j^^, etc. denote mean values of the 
corresponding operators. 
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mean EM field is not weak. So, we have to separate the condensate mode from 
the photon degrees of freedom in the Hamiltonian and the statistical operator. 

4.1 The time- dependent unitary transformation 

The condensate mode is most easily isolated by introducing the r-dependent 
unitary transformation 

g{n, t) = e*^("'") g{n, r) e-*^("'^), (4.1) 

where the operator C {n, r) is given by 

C{n,T) = jda {A'i{x)Ii^^{xJ-Ii,^{x)A^A^,)]. (4.2) 

Note that the unitary transformation (4.1) does not affect fermionic operators 
and has the properties 

giC(n,r) ^^(3.^) Q-iC{n,T) _ _^ ^/.(^^j^ 

giC(n,r) Q-iC{n,r) ^ ^jm^^^^, ^ jj^^x). (4.3) 

Taking now into account that, for any operator O, 

{Oy = Tr {e^^("'^) 6 e-*^("'") ^(n, r)} = (e*^^"'^) 6 e'^^^"'^))^ , (4.4) 

we find that 

(i:^(a;j)J = (n^^(a;j)J = 0. (4.5) 

Thus, in the state described by the transformed statistical operator (4.1), the 
canonical dynamical variables and fix have zero mean values and, hence, 
they are not related to the condensate mode. In other words, after the unitary 
transformation the EM field operators correspond to the photon degrees of 
freedom. Based on the above arguments, it is convenient to use ^(n, r) as the 
statistical operator of the system. It should be noted, however, that Q{n. r) 
does not satisfy the von Neumann equation (2.8) since the operator C depends 
on T. In order to derive the equation of motion for g{n,T), we differentiate 
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Eq. (4.1) with respect to r. After some algebra which we omit, we find that the 
transformed statistical operator satisfies the modified von Neumann equation 



dg{n, ■ 



— I 



g{n,T),n^{n) 







(4.6) 



with the effective Hamiltonian 



7Y^(n) = e^^("'^) H^{n) e"^^^"'^) 




^i«.j-MMn.>j|. (4.) 



The transformation of H{n) in the first term is trivial due to Eqs. (4.3) and 
the fact that the transformation does not affect fermionic operators. It is 
convenient to eliminate the derivatives in the last term of Eq. (4.7) with the 
aid of the equations of motion for the condensate mode 

^ = -n'^(a:), ^ = V,F^(x) - ,-(.), (4.8) 



which are easily derived using Eqs. (2.23), (2.24), and the canonical commu- 
tation relations (3.21). The tensor in Eq. (4.8) is the mean value of the 
operator (3.37), and j±{x) is the transverse part of the mean polarization 
current 

fix) = {j\x,)Y. (4.9) 

Inserting Eqs. (4.8) into Eq. (4.7), the effective Hamiltonian can be written 
as a sum 

7i»=7YS(n)+7irntW- (4-10) 

The main term 

niin) = Hj,{n) + Hem + I daj^{xj A^{x) (4.11) 

describes free photons and fermions interacting with the total electromagnetic 
field 

A''{x)=AUx)+A^{x), (4.12) 
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where the mean polarization field A'^{x) is given by 



) 



T 



(4.13) 



The term Hl^^ln) in Eq. (4.10) describes a weak interaction between photons 
and fermions. The explicit expression for this term is 



+ 2/^^/ rf(T'Ai||(a;,;r)G(a;,-a;'JAi|,«;r), (4.14) 



where the operators 



(4.15) 



represent quantum fluctuations of the fermionic current. The essential point 
is that now the interaction term (4.14) does not contain a contribution from 
the condensate mode and, consequently, one can use perturbation expansions 
in the fine structure constant. 

4.2 Maxwell equations 

To complete our discussion of the condensate mode, we will show how the 
Maxwell equations for the total mean EM are derived in our approach. Ac- 
cording to Eq. (4.12), the total field tensor can be written as 



The external field tensor F^^^ is assumed to satisfy the Maxwell equations 



with some prescribed external current j^xt- On the other hand, the polarization 
field tensor F^'^[x) is the mean value of the operator 



(4.16) 



5,i^e'it(^)=j:xt(^) 



(4.17) 



F^'''{x^)^^^'A'' -d^'A" 

= Fr + rf {pi + VM||) - {p + Vi,,) . 



(4.18) 
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Recalling Eqs. (2.23) and (2.24), straightforward algebraic manipulations with 
the equations of motion for the field operators lead to the Maxwell equations 
for the polarization field tensor 



Now Eqs. (4.17) and (4.19) can be combined into the Maxwell equations for 
the total field tensor 



A solution of these equations gives the total mean field A'^ in terms of the 
total mean current. 



5 Kinetic description of QED plasmas 

5.1 The "one-time" Wigner function 

Within the hypcrplane formalism a natural way of describing kinetic processes 
in the fermion subsystem is by the "one-time" Wigner function which depends 
on the variable r. Since there is the gauge freedom for the mean field A'^, it is 
convenient to employ the gauge-invariant Wigner function on the hyperplane 
(T„ defined as 



(4.19) 



d^J^^''{x)=f{x)+jUx). 



(4.20) 




X exp <ieA{Xj^ + 



b.-^r) (5.1) 



with the gauge function 



X 




1 




(5.2) 







In Eq. (5.1) the one-particle density matrix is the mean value 



Paa' ix^,x'^;r) = ix^,x'jy = ix^,x'jy- 



(5.3) 
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of some density operator p^,. In the literature one can find different definitions 
for the fermionic density operator. The most often used definitions are 

Paa'(3^x>0 = -^[V'a(^x)>^a'«)]> (5-4) 
Pla'(2^x>0 ^^^a'lOV'aW: " (5-5) 

These two operators are related by 

P<m'{x^, X'j) = p'aa'ix±, x'^) + ^aa'(^x, ^'J, (5-6) 

where the last c-number term represents the vacuum expectation value of p 
since the vacuum expectation value of p' is zero. It can be shown, however, 
that the vacuum term in Eq. (5.6) does not contribute to local observables 
like the mean current j^{x). The advantage of the definition (5.5) is that the 
mean values of one-particle dynamical variables (summation over repeated 
spinor indices) 

6 ^ j dada'0^,^{x'^,xj ■'ta'M '^a{x±) ■ (5-7) 



are conveniently expressed in terms of the density matrix p' — {p'Y: 

{Oy = j dada' 0^,^{x'^,xJp'^^,{x^,x'^;t). (5.8) 



Unfortunately, the equation of motion for the density operator (5.5) with the 
Hamiltonian (4.10) contains vacuum (divergent) terms. On the other hand, 
such terms do not appear in the equation of motion for the operator (5.4). 
For this reason, we shall take the operator (5.4) as the one-particle density 
operator in Eq. (5.3). An analogous definition was used previously for the 
phase-space description of the QED vacuum in a strong field [7] . 

The Wigner function (5.1) is defined on a given family of hyperplanes cr„ ^ and, 
hence, depends parametrically on the normal four- vector n. It should be noted, 
however, that local observables calculated from the Wigner function do not 
depend on the choice of n. As an important example, the mean polarization 
current (4.9) can be written in the form 



f{x)^e{:i^ixjri;{xj:y 

= e / -0^S{p-n)tr{rW{x^,p^;r^x-n)}, (5.9) 
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where the symbol "tr" stands for the trace over spinor indices. Geometrically, 
the above relation means that, in calculating the current, the invariant time 
T has a value such that the space-time point x lies on the hyperplane cr„^. 



5.2 The photon density matrix 



To define the photon density matrix, we start from the plane wave expansion 

of the vector potential operator in terms of creation and annihilation op- 
erators. By analogy with the well-known representation for the free photon 
field in the special Lorentz frame where = (1, 0, 0, 0), we write 



X E {«;(?x) e-^''-^ + a\{q.) e^''"} , (5.10) 

1=1,2 

where e^{q^,l) are real- valued polarization four- vectors and 

is the dispersion relation for free photons on the hyperplane. The conditions 
V^^x = nij^A'i — mean that the polarization vectors satisfy 

q.,e>^iQ.,l)^n^e'^{q,,l)^0. (5.12) 



The expansion of the operator fix into plane waves is found from (5.10) by 
using = -A'i: 



J -/2a;,(gJ(27r)3 ^ ^ 



X E e'^(?x, e-^^-^ - ajiqj e^^-} . (5.13) 



1=1,2 



Assuming the commutation relations for the creation and annihilation opera- 
tors 



(5.14) 

K(?x),%(5x)] = K(9x),4(9x)] = 0, 
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and the completeness relation for the polarization vectors 

E e^(?x, e^(?x, = - i^'" - ^-jf] , (5.15) 

1=1,2 \ ^± / 

the commutation relation (3.21) for the field operators is recovered. Note that 
Eqs. (5.10) and (5.13) give the field operators in the interaction picture. The 
corresponding expansions for the field operators in the Schrodinger picture are 
obtained by setting r = 0. In this case the delta-function S{q^^ — uj^{q^)) can 
be replaced by S{q^^). 

The above considerations suggest that it is natural to define the photon density 
matrix in terms of the Schrodinger operators 

which satisfy the commutation relations 

The photon density matrix is defined as 
where 

Nii'{x^, x'J = (fiUx'J 0i{xj (5.19) 

is the photon density operator. It should be emphasized that in Eq. (5.18) 
the average is calculated with the transformed statistical operator g{n, r) in 
which the condensate mode of EM field has been eliminated. When written in 
terms of the average with the statistical operator g{n, r) , the photon density 
matrix takes the form 

NAx,,x',;t) = (4'(^x,0)' - {0i{xs.)y (4(^'x)>^ (5-20) 
where the last term corresponds to the contribution from the condensate mode. 
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^(g„)e-«-a,(gJ, 



/ (27r)3/2 
/7^%„)-"-^«K^x), 



(5.16) 



(5.17) 



(5.18) 



5.3 The covariant statistical operator in QED kinetics 



The evolution of the fermionic Wigner function (5.1) and the photon density 
matrix (5.18) is governed by kinetic equations which can be derived from the 
equations of motions 



^ Paa' (^x, r) = -i lY { [/)„„, {x,, x'J,ni{n) + K[,,{n)] ^(n, r) } , (5.21) 

|:A^,,(a;,,<;r) = -iIv{[iV,,(x,,a;l),7i;S(n)+K„,(n)]^(n,r)}. (5.22) 

There are two ways to express the right-hand sides of these equations in terms 
of the fermionic and photon density matrices using perturbation expansions 
in the fine structure constant. One method is by considering the hierarchy 
for correlation functions which appear through the commutators with the 
interaction Hamiltonian ?Y^t(n) and then employing some truncation proce- 
dure. Another method is to construct an approximate solution of Eq. (4.6) 
in terms of the density matrices p and A^. In both cases one has to impose 
some boundary conditions of the retarded type on the correlation functions 
or the statistical operator. The standard boundary condition in kinetic theory 
is Bogoliubov's boundary condition of weakening of initial correlations which 
implies the uncoupling of all correlation functions to one-particle density ma- 
trices in the distant past, i.e., for r — oo. In the scheme developed by 
Zubarev (see, e.g., [21]), such boundary conditions can be included by using 
instead of Eq. (4.6) the equation with an infinitesimally small source term 



dQ{n, t) _ . 
dr ' 



Q{n,T),lrr{n) = -£ {^(n, r) - ^,gi(n, r)} , (5.23) 



where £ — > -|-0 after the calculation of averages. Here Q^^-^in^^T) is the so- 
called relevant statistical operator which describes a Gibbs state for some given 
nonequihbrium state variables. In QED kinetics these variables are the Wigner 
function (5.1) and the photon density matrix (5.18). Therefore, following the 
standard procedure [21], we obtain the relevant statistical operator in the form 
(with summation over spinor and polarization indices) 



Qr^\in,T) = — exp<j - / dada' Ai2(a;^, x'^; r) :'0„(a:x) V'a'Ca^l) : 




(5.24) 
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where Z^^i(n, r) is the normalization constant (or the partition function in the 
relevant ensemble) and X[fJ,{x^, x'^; r), \\i,''\x^,x'^; r) are Lagrange multipliers 
which are determined by the self-consistency conditions 

(5.25) 

Nii'{x^,x'^; r) = IV {Nii,{x^, x'Jg,^i{n, r)} . 
Using Eq. (5.23) for the transformed statistical operator leads to the hierarchy 



- £ {^i-fe(3^ix, • • • , ^fcx; ^) - {Oi{.x^^) ■ ■ ■ Ok{xkS)\ J , (5.26) 

where 0^{x^^ are some Schrodinger operators which may depend on the 
fermion operators as well as on the EM operators, and 

K-ki^i^^ Xk±; n, t) = (Oi(xi^) • • • Ofe(Xfe^))^ (5.27) 



are the "equal-time" correlation functions in which the condensate mode of 
the EM field is eliminated. Since the relevant statistical operator (5.24) ad- 
mits Wick's decomposition, the last term in Eq. (5.26) ensures the boundary 
condition of complete weakening of initial correlations. Note that the explicit 
knowledge of the statistical operator ^(n, r) is not needed when considering 
the hierarchy for the correlation functions. Use of some truncation procedure 
is a standard practice in this case. The hierarchy for correlation functions will 
be discussed in subsequent papers in the context of the derivation of collision 
integrals. 

Another method of handling Eq. (5.23) is by considering its formal solution 

T 

Q{n, r)^e I dr' e'^^^"^') U(r, r') g,,,(n, r') U\r, r'), (5.28) 

—oo 

where the evolution operator can be written as the ordered exponent 

U{t,t') = exp I -i Jn^{n) df\ . (5.29) 
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After partial integration, the expression (5.28) becomes 



g(n, r) = g^^i{n, r) + A^(n, r), 



(5.30) 



where 



T 



Ag{n,T) 




— oo 




dr' 



i[g,,,{n,T'),n^'{n)]\u\T,r'). (5.31) 



The representation (5.30) for the statistical operator allows to separate the 
mean-field terms and the collision terms in Eqs. (5.21) and (5.22). Taking into 
account the self-consistency conditions (5.25) and the fact that the Hamilto- 
nian (4.11) is bilinear in the fermion and photon operators, we arrive at the 
equations 



|-iV,,(x-„<;r) = -t{[Nu,ix,,x'J,HEj^]y + lif'\x,, x',;t), (5.33) 



In the presence of a strong EM field, the evolution of the fermion subsystem 
is governed predominantly by its interaction with the mean EM field. Thus, 
the covariant mean- field kinetic equation for the Wigner function (5.1) can 
be derived from Eq. (5.32) neglecting the collision integral. This kinetic equa- 
tion as well as the collision integrals (5.34) and (5.35) will be considered in 
subsequent papers. 



where the collision integrals for fermions and photons are given by 



/^''\x„x1;t) = -tTr{[Ni,{x,,x'J,nUn)]Ag{n,T)}. 



(5.34) 



(5.35) 
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6 Concluding remcirks 



We have shown that the hyperplane formahsm can serve as the basis for kinetic 
theory of QED plasmas in the presence of a strong external field. The formal- 
ism has the advantage that it is manifestly covariant and therefore allows 
to introduce different approximations in covariant form. Only minor changes 
with respect to the non-relativistic density matrix method are introduced, so 
that many well-developed approaches can be directly applied to QED plas- 
mas. For instance, the explicit construction of the statistical operator allows 
to incorporate many-particle correlations through the extension of the set of 
basic state parameters (see, e.g., [21]). Note also that, using the Heisenberg 
picture on hyperplanes, nonequilibrium Green's functions can be introduced 
with respect to the invariant time parameter r. In such a way, the spectral 
properties of microscopic dynamics can be incorporated. 

The scheme outlined in this paper is also applicable to other field theories, 
like QCD transport theory. In QCD, however, additional problems arise due 
to its non-Abelian structure, which needs further considerations. 

Finally, we would like to emphasize once again two key problems in a covariant 
density matrix approach to relativistic kinetic theory in the presence of a 
strong mean field. First, it is necessary to perform canonical quantization of 

the system on a hyperplane in Minkowski space. Second, the condensate mode 
must be separated from the quantum degrees of freedom at any time. We have 
shown how these problems can be solved in the context of QED plasmas. As a 
result, a general form of kinetic equations for fermions and photons was given. 

The scheme outlined in this paper is also applicable to some quantum field 
models used in QCD transport theory. In this case the non-Abelian algebra 
must be worked out to describe the quark-gluon plasma. 



Appendix A 

Commutation relations for electromagnetic field on hyperplanes 

Let us write the constraint equations for the canonical variables A± and 11' 
in the form Xn{^±) = 0) where 



(A.l) 
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For any functionals and $2 '^^ the field variables Aj_ and Uj^, we define the 
Poisson bracket 

where the constraints are ignored in calculating the functional derivatives. 
Applying this formula to the canonical variables we obtain 

[^^^(xj,n,,«)]p = 6'^j\x,-x'J (A.3) 

with the three-dimensional delta function (3.16). All other Poisson brackets 
for the canonical variables are equal to zero. In the Dirac terminology, func- 
tions (A.l) correspond to second class constraints since the matrix 

is non-singular. A straightforward calculation of the Poisson brackets shows 
that the non-zero elements of C are 

C,^{x^,x'J = -C^i{x^,x'J = -V^V^S\x^ - x'J, 

C^^ix^, x'J = -C^six^, x'J = 5^ix^ - x'J. (A.5) 

According to the general quantization scheme [18,19], commutation relations 
for canonical operators are defined by the Dirac brackets for classical canonical 
variables. In our case the Dirac brackets are written as 

- / o^f^ / d(^'[^i,XN(x±)]pCNN'(x^,x'J[xN'W±),^2\p (A.6) 

(summation over repeated indices). The inverse matrix, C]^lf,{x^, x'^^, satisfies 
the equation 

J d(T C jqj^ii{x ^,x ^ C j\fii]\^i{x ^,x ^ — 5jqj^i 5^{x ^ ^x)- (-^■'^) 

Since the matrix elements (A.5) of C depend on the difference x^ — x'^, 
Eq. (A. 7) can be solved for C^^ using a Fourier transform on cr^^, which 
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is defined for any function f{x) as 

/(r,pj = J d^xe'P-''6{x-n-T)f{x). (A.8) 



The inverse transform is 

fix) ^ f{T,xJ = I ^,e-^'^-^S{p-n)f{T,pJ. (A.9) 



If we perform the Fourier transformation in Eq. (A. 7), we find by insert- 
ing (A. 5) that the non-zero elements of C"^ are 



C^;{x,,x'J = -C^,\x,,x'J = -j ^"^"■^'"''^^(^^•^)^' 
C3-/(x,,x'J = -C4-/(x„x'J = -5\x,-x':). (A.IO) 

Now the Dirac brackets (A. 6) for the canonical variables are easily calculated 
and we obtain 



K(a;J,n::«)]j, = c'^^^x^-x'J, (A.ll) 

aKx'JIj, = [n'^(xj,ni(a;'j]i, = o, (A.12) 

where the functions d^^{xj^ — x'^ are given by Eq. (3.23). According to the 
general quantization rules, the commutation relations for canonical operators 
correspond to i[. . .j^. Thus, in the hyperplane formalism, the commutation 
relations for the operators of EM field are given by (3.21) and (3.22). Obviously 
these relations are valid in the Schrodinger and Heisenberg pictures. 



Appendix B 

Anticommutation relations for the Dirac field on hyperplanes 



To find the anticommutation relations for the fermion operators on the hyper- 
plane cr„,^, it is sufficient to consider a free Dirac field. Our starting point 
is the standard quantization scheme in the frame where x'^ = (t,r) and 

— (1,0,0,0) (see, e.g., [17]). In that case the field operators ip^ and ip^ 
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can be written in terms of creation and annihilation operators according to 

/• d'^p 5{jP — e{p)) 



where e{p) = yp'^ + m? is the free fermion dispersion relation. Constructing 

the expression ^/^^/(a;')} for two arbitrary space-time points and recall- 

ing the anticommutation relations 

{hs{p)Mp')} = {4(P),4(P')} = 5ss'5\p-p'l (B.l) 



as well as polarization sums 

E «as(p)«a'sb) = + m 



s=±l 

we arrive at 



s=±l 



7% - m 



V..W.4M}=/g^{K + 



m 



-ip-{x—x') 



+ [7% 



m 



J,p-{x—x 



where p'^ — yp^ + m?- Using 
(fp 1 



(27r)3 2e(p) J (27r) 



pO>0 



(B.2) 



(B.3) 



Eq. (B.2) can be rewritten in a Lorentz invariant form 



{V^„(x),|„,(x')} = / ^ {[7'^p^ + m]^^,e-^-(---')5(/-m^ 



(27r)= 



pO>0 



+ 



TP^ - m 



Jpix-x') r/ 2 



S{p — m 



pO>0 



(B.4) 



The anticommutation relation on the hyperplane ^ is now obtained by set- 
ting X — nT-\-x_^ and x' — nT-\-x'_^. In calculating the integrals, it is convenient 
to use the decomposition = n'^p,, -\-p'i, (p., > 0). Then we get 
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X { [7||P|| + 7>x^ + m]^^, e-ip^,i<-^'D 



+ 



7|,P„ + llP., - m] ^, } (B.5) 



with the dispersion relation on the hypcrplanc 



e(Px) = sj-V.^^. + VP?. (B.6) 

Finally, changing the variable px —p± in the second integral in Eq. (B.5), 
we obtain the anticommutation relation (3.24). The relations (3.25) can be 
derived by the same procedure. 
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